Experimental Evidence for Two-Dimensional Magnetic Order in Proton Bombarded 

Graphite 
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We have prepared magnetic graphite samples bombarded by protons at low temperatures and 
low fluences to attenuate the large thermal annealing produced during irradiation. An overall 
optimization of sample handling allowed us to find Curie temperatures T c > 350 K at the used 
fluences. The magnetization versus temperature shows unequivocally a linear dependence, which 
can be interpreted as due to excitations of spin waves in a two dimensional Heisenberg model with 
a weak uniaxial anisotropy. 
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Recent advances to develop nanographitic systems 
have led to a renewed interest on their electrical prop- 
erties worldwide [H- A single layer of graphite, the two- 
dimensional (2D) graphene, appears to have quantum 
properties at room temperature Q as well as rectifying 
electronic properties [H, Q|. On the other hand, some 
of those properties were already observed in highly ori- 
ented pyrolytic graphite (HOPG) of low mosaicity, as the 
quantum Hall effect H and de Haas - van Halphen quan- 
tum oscillations even at room temperature Q. The two- 
dimensional properties of the graphene planes in graphite 
open up the possibility of using nanometer to micron size 
regions of graphite in new integrated devices with spin- 
tronic properties either through the use of ferromagnetic 
electrodes, e.g. spin-valves, and/or making graphite it- 
self magnetic. In fact this has been a topic of study in the 
last years and reports exist showing magnetic hysteresis 
in blank graphite but especially in proton bombarded 
graphite [8J. Severe limitations in the sensitivity and re- 
producibility of standard magnetometers added to an- 
nealing effects during bombardment, hindered the identi- 
fication of a critical temperature T c as well as the charac- 
teristics and dimensionality of the ferromagnetic signals. 
The aim of this work is to show that specially prepared 
highly oriented pyrolytic graphite (HOPG) samples show 
ferromagnetic order with T c > 350 K and the magnetiza- 
tion temperature dependence is in good agreement with a 
2D anisotropic Heisenberg model (2DHM) and the pres- 
ence of spin waves excitations 0, 0, EH • 

For the experiments we used four pieces of a HOPG 
sample grade ZYA, samples 1 to 4 (mass: 12.8, 12.5, 
10.1, and 6 mg respectively) irradiated by a 2.25 MeV 
proton micro-beam (sample 4: 2.0 MeV, 0.8 mm broad 
beam) perpendicular to the graphite planes. With the 
micro-beam we produced several thousands of spots of 
~ 2 /im diameter each and separated by 5 /im (sam- 
ple 1) or 10 /im (samples 2 and 3) distance, similarly to 
the procedure used in Ref. 12. Samples 1 and 2 were 
irradiated at 110 K whereas samples 3 and 4 at room 
temperature. Further irradiation parameters for sam- 



ple 1 (2,3,4) were: 51375 (25600,25600,6) spots, fluence: 
0.124 (0.08,0.13,0.3) nC//im 2 , total irradiated charge 
46.9 (44.8,37.4,900) fiC, and 1 nA proton current (100 nA 
for sample 4). The pieces we have irradiated showed an 
iron concentration (the only detected magnetic impurity) 
within the first 35 /xm of ~ (0.4±0.04) fig/g (< 0.1 ppm). 

Previous experiments [8( showed ferromagnetic mag- 
netic moments at saturation m sa t ~ 1 /iemu and there- 
fore put severe constrains to experimentalists, not only 
regarding the sensitivity of the used magnetometer but 
also its reproducibility after sample handling. In this 
work two main experimental improvements have been 
achieved. Firstly, we enhanced the ferromagnetic part 
produced by irradiation reducing annealing effects. In 
samples 1 and 2 the micrometer spots were produced 
at a nominal temperature of 110 K during irradiation 
(18 hours). For comparison and to reduce further anneal- 
ing effects sample 4 was irradiated with a broad beam 
and low fluence. Second, we have designed a sample 
holder that allows us to measure the magnetic moment 
of the sample in the SQUID and to fix it inside the irra- 
diation chamber without any changes. We investigated 
the reproducibility of the magnetic measurements and 
checked that the sample holder handling (with sample) 
[l3j |. i.e. inserting it and taking it out of the irradia- 
tion and SQUID chambers 14|, does not produce sys- 
tematic changes of the magnetic signal. Our arrange- 
ment provides a reproducibility of ~ 10~ 7 emu in the 
measured field range and allows the subtraction of the 
virgin data from those after irradiation point by point, 
increasing substantially the sensitivity of the magnetic 
measurements to ~ 2 x 10~ 8 emu. 

Figure 1 shows the hysteresis loops of the magnetic 
moment m of sample 2 at two temperatures. These loops 
are obtained directly from the difference of the measure- 
ment "after" minus "before" irradiation. The loop at 
5 K as well as the measured temperature dependence 
at constant field indicate a paramagnetic contribution 
nip = 0.575 H/T[/j,eimi K/kOe] for this sample, i.e. less 
than 10% of the ferromagnetic signal at 3 kOe. At 300 K, 
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FIG. 1: Magnetic moment as a function of applied field for 
the irradiated sample 2 at 300 K (o) and 5 K (■) obtained 
after subtracting the data of the non irradiated sample. The 
points (*) are obtained for the same sample at 300 K after 
taken out the first ~ 5 micrometers from the irradiated surface 
side. The inset blows out the data at low fields to show the 
finite hysteresis and the clear temperature dependence of the 
coercive field and remanent magnetic moment. 



however, m p is negligible. These loops, their tempera- 
ture dependence as well as the finite hysteresis, see inset 
in Fig. 1, indicate the existence of magnetic order with 
a Curie temperature higher than room temperature. 

Sample 3, which was irradiated with similar number 
of spots, fluence and total charge but at room temper- 
ature, shows a a ferromagnetic signal at saturation ~ 5 
smaller than that obtained for samples 1 or 2, in agree- 
ment with previous work|8[. These results indicate the 
reliability and sensitivity of the used procedure as well 
as the absence of obvious artifacts in the measurements. 

After peeling out the first micrometers from the irradi- 
ated surface of sample 2 the ferromagnetic contribution 
decreased by one order of magnitude, see Fig. 1. We can 
answer now the question whether the Fe concentration 
in the sample and due to some hypothetical annealing by 
the protons could be responsible for the observed ferro- 
magnetic signal. In the first micrometer and taking an 
irradiated area < 0.026 cm 2 , the magnetization at room 
temperature is then > 0.5 emu/g. In this region we esti- 
mate that the mass of the ferromagnetic carbon material 
is < 6 fig. Were the measured Fe concentration ferromag- 
netic at 300 K then it would contribute with a magnetic 
moment < 0.6 x 10~ 10 emu, i.e. 50.000 times smaller 
than the measured one. With the mass of the ferromag- 
netic part of the irradiated HOPG sample we estimate a 
magnetic moment per carbon atom mc ^ 0. 001 (is, in 
very good agreement with XMCD results 15 1. 

Figures 2 and 3 show the temperature dependence of 
the ferromagnetic moment for samples 1 and 4, respec- 
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FIG. 2: Normalized magnetic moment (m(0) = 2.60 /iemu) 
obtained for irradiated sample 1 at 2 kOe. The data points 
are obtained after subtracting the data from the sample before 
irradiation and a paramagnetic (Curie) contribution m p (T) = 
4.9/T[/iemu K]. The errors bars indicate typical errors due 
to the subtraction of the data from the virgin sample. The 
chosen parameters for the theoretical curves are T c — 360 K, 
T° w = 850 K (A = 0.001). The continuous line is obtained 
from (JTJ|. The dotted line is the 3D Bloch T 3/2 model with 
spin waves [3]. The dashed-dotted line with close triangles 
shows the results of a Montecarlo simulation with anisotropy 
(square lattice of 200 x 200 points). 



tively. Because the paramagnetic signal contributes sig- 
nificantly only at T < 25 K, we have subtracted it in 
both figures in order to show only the ferromagnetic part. 
Up to the highest measured temperature of 380 K this 
magnetic moment behaves reversible. Furthermore, no 
changes in m within experimental error were observed 
after leaving the samples several months at room tem- 
perature. 

One of the interesting and indicative results shown in 
Figs. 2 and 3 is the unequivocal linear dependence. This 
is an indication of 2D magnetism and the slope can be 
interpreted as due to the excitation of 2D sgin waves 
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that reduce the magnetization linearly with T 
We are not aware of any model hamiltonian producing 
such a linear behavior in m(T). Therefore, to analyze the 
measured temperature dependence we discuss the 2DHM 
with anisotropy that provides a linear dependence with 
T. The discrete Hamiltonian describing the 2DHM reads 
H = -JJ2iji s iz S ]z + (1 _ A)(Si X S JX + SiySjy)}, where 
Si = (Si X , Sty, Siz) represents a unit vector in the direc- 
tion of the classical magnetic moment placed at the site i 
of a 2D lattice. The sum is performed over all near- 
est neighbor pairs, and J is the exchange coupling. The 
parameter A represents the uniaxial anisotropy in the z- 
direction. The case A = is the isotropic 2DHM and is 
known to have T c — 0. However, just a small anisotropy 
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raises T c considerably because T c ~ — 1/ In A for A — > 0. 
It can be shown [ij, [l^! 11 1 that the normalized spin- 



waves magnetization in the anisotropic axis behaves as 
M s z w = 1 - T/T* w - 2T 2 /(T*T* W ) - (2/3) (T/Tf) 3 at 
low temperatures, where T* = 4J. This result is ob- 
tained using perturbation theory techniques [3, [13] up 
to third order in spin waves. The parameter T^ w is the 
spin wave critical temperature due to low-energy spin 
wave excitations; it is given by k B T^ w = 2n J/K(l — A), 
where K(x) is the elliptic function. Near the critical 
temperature T c the physics can be better described by 
a 2D Ising model that should provide a good descrip- 
tion of the spin flip excitations. Then T c is given by 
T C (J) = 1.52 where J is the renormalized exchange 
due to the spin waves excitations according to the ex- 
pression J(T) = J[l - 2T/T c stu ]. The values of M z at 
T < T c can be expressed as: 



M Z {T) « Mf (T, J)M*[T, J(T)] 



(1) 



The first factor in the rhs of |TJ) is the magnetization 
due to spin waves and the second one is the magneti- 
zation due to an Ising model with the exchange renor- 
malized by the spin waves. We have checked this theo- 
retical result against Montecarlo calculations with A = 
0.001 and the agreement is excellent, especially at low 
anisotropies ll| as it is shown in Figs. 2 and 3. In Fig. 2 
we have plotted also the normalized spin waves contribu- 
tion M* w /M* w (0) up to third order. The Heisenberg re- 
sult approximated by |T]) and the Montecarlo calculation 
agree and both fit the experimental data with the param- 
eters T° w = 850 K, J{T C = 360K) = 237 K, indicating 
an anisotropy A ~ 0.001. Sample 2 shows a similar be- 
havior and its data can be fitted with T" w ~ 1000 K, 
J(T C ~ 310K) = 202 K. The data for sample 4 shown 
in Fig. 3 show also a linear behavior. Extrapolating 
the SW contribution to m(T*) ~ we conclude that 
T c < T* ~ 640 K. Then using (1) we estimate T c > 450 K 
with A < 10~ 4 , see Fig. 3. These results already show 
that T c increases with fluence, provided that one can re- 
duce simultaneously the annealing effects produced dur- 
ing irradiation. For comparison we also have plotted in 
Figs. 2 and 3 the Ising model result that has no spin 
waves and the 3D Bloch T 3 / 2 law that includes spin waves 
[l8[ . The comparison indicates clearly that spin waves in 
2D dominate the magnetization up to > 300 K and that 
the usual 3D model does not fit the data. 

There is no doubt that defects in the graphite struc- 
ture are one of the possible origins for localized magnetic 
moments. The ferromagnetism triggered by the bom- 
bardment should be correlated to the produced defects 
located at approximately the first micrometer from the 
sample surface. To discuss a mechanism responsible for 
the coupling between the magnetic moments, we need 
first to estimate the density of defects. For sample 1 we 
have 0.9 nC total irradiated charge per spot in an area 
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FIG. 3: Normalized magnetic moment (m(0) — 4.9 /iemu at 
10 kOe) obtained for sample 4 at (10,3,1) kOe (•,□,★) after 
subtracting the data from the sample before irradiation and a 
paramagnetic (Curie) contribution m p (T) = 1.18H /T[fj,emu 
K/kOe]. Different theoretical curves are the same as in Fig. 2 
but with parameters: T c = 450 K and T* w = 1050 K. 



of (~ 7r0.6 2 )/Ltm 2 . Using SRIM2003 Monte Carlo simula- 
tions with full damage cascades and 35 eV displacement 
energy we obtain a vacancy density of ~ 5 x 10 20 cm' 3 
at the surface, which means a distance between vacancies 
I ~ 1.3 nm~ 9a, where a = 0.14 nm. This distance is 
much smaller than the inverse of the Fermi wave vector 
1/fcp ~ 30 nm for a Fermi energy of 20 meV or calculated 
using the 2D carrier density 0]. 

Regarding the coupling needed to have room temper- 
ature magnetic ordering there is in first place the di- 
rect coupling for nearly localized spins at the defects, 
which should be in the range of ~ 300 K. Recently the 
RKKY coupling between large defects in graphene has 
been studied for Fermi energy tending to zero[l9|. This 
coupling might be always ferromagnetic because fci?r <C 1 
for r ~ I. However, estimations of the Curie tempera- 
ture for this coupling within our defect densities provide 
values of the order of 20 K. What appears important 
is a super-exchange mediated by the two different sites 
in the graphite lattice 0, E} or between magnetic mo- 
ments from defects and from hydrogen atoms, which may 
effectively increase the magnetic moment density on a 
graphene lattice. 

We note that large concentration of hydrogen is found 
in the first micrometer thick region at the surface of 
graphite samples j22|. Therefore we should take into ac- 
count the possible influence of hydrogen in triggering lo- 
calized as well as non-localized magnetic moments in the 
graphite layers 2(J 23| . Irradiation may contribute as de- 
fect generation as well as dissociating the existing molec- 
ular hydrogen enabling its diffusion and bonding in de- 
fective parts of the lattice structure. All these moments 
will tend to be ferromagnetically coupled enhancing the 
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Curie temperature by the RKKY coupling. 

Within this picture it becomes clear that the enhance- 
ment of the defect density, which occurs at larger depths 
from the surface in the inner part of the irradiation path 
up to full amorphization at a depth ~ 35 ... 40 /im, 
perturbs too much the graphene lattice destroying in 
this way the necessary band structure and carrier den- 
sity. This may explain the experimental observation of 
a rather well-defined critical temperature (and not a dis- 
tribution) and also the difficulty one has to reach much 
higher ferromagnetic magnetization values increasing the 
proton fluences clearly above the values used here. If an 
electron-mediated coupling between defects plays a role, 
we expect that for an adequate defect density it should 
be possible to influence the magnetic order shifting the 
Fermi energy by applying an appropriate bias voltage. 

The results of samples 1 and 2 provide clear evidence 
for the good reproducibility of our approach: although 
the spot density, beam diameter as well as total charges 
were different, the produced defect densities in the irra- 
diated paths were similar for both samples and therefore 
we expect to obtain similar critical temperatures as the 
measurements showed. Changing the defect density as 
well as their distribution in the lattice one may tune the 
ferromagnetic transition temperature as well as the mag- 
nitude of the magnetization produced by irradiation, as 
the data for sample 4 clearly indicate. As a rule of tumb 
robust ferromagnetism with T c > 300 K by proton irra- 
diation in graphite can be reached with fluences of the 
order of 0.1 nC//im 2 . 

In conclusion our work shows that irradiation of mi- 
crometer spots in graphite at low temperatures as well 
as broad irradiation, both at very low fluences, increases 
significantly the magnitude of the magnetic order with 
Curie temperatures T c > 300 K. The use of especial sam- 
ple holders made possible to reduce sample handling be- 
tween irradiation chambers and SQUID measurements to 
a minimum, ruling out simple introduction of impurities 
or the influence of operative artifacts. This approach in- 
creased substantially the sensitivity and reproducibility 
of the magnetization measurements allowing us to ob- 
tain directly the effects produced by irradiation within 
an error of ~ 10~ 7 emu. The experimental localization 
of the ferromagnetic irradiated part of the sample indi- 
cates that the graphite structure is important and that 
at the used proton energies low fluences are preferential 
to trigger a robust ferromagnetic order. We showed that 
the magnetization of the magnetically ordered contribu- 
tion decreases linearly at T < T c , a behavior that can 
be assigned to the signature of low energy spin waves ex- 



citations well described by an uniaxial two dimensional 
anisotropic Heisenberg model. 
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